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Abstract: The temperature dependent BCS equations are modified in order to 
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I. INTRODUCTION 

in : 

The effect of temperature upon pairing correlations in nuclei was studied long ago [|TJ, 
CN ' after the theory of superfluidity was introduced to describe the appearance of a gap in the 
low-energy spectrum of finite nuclei 0. The extension of the theory of superfluidity to the 
Q\ • finite temperature case was prompted by the search of evidence about the nuclear phase 
transition from a superfluid to a normal state, in analogy with the one found in condensed 
\ matter systems ||. The vanishing of pairing correlations with temperature is related to the 

fact that the excitation energy breaks pairs of particles which block the single particle levels 
close to the Fermi surface, where the pairing correlations are initiated. 

The change of the pairing properties with the temperature was also studied in competi- 
tion with the angular momentum In this case a new effect was observed, which is due 
to the fact that at zero temperature the single-particle states close to the Fermi level are al- 
ready partially blocked by the un-paired nucleons which form the finite angular momentum 
of the nucleus. Thus, when the temperature is switched on the first effect is a depletion of 
the partially blocked single-particle states which, in turn, induces an enhancement of the 
pairing correlations. This behaviour is different for the case of thermal excitations of nuclei 
with zero angular momentum, where the pairing correlations decrease monotonically with 
the increase of the temperature. 

The temperature effects on pairing correlations were studied both in the BCS |T|^J^,0] and 
HFB || approximations, usually with a constant single-particle level density approximation. 
The available theoretical evidence shows that for medium and heavy-mass nuclei the pairing 
correlations disappear for a critical temperature of the order of 0.5-1.0 MeV. The excitation 
energies corresponding to this critical temperature are quite small for nuclei close to the 
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beta stability line and therefore in all such calculations the coupling with the continuum 
spectrum was neglected. The situation is different for nuclei which are far from the stability 
line. In this case the Fermi level lies close to the continuum threshold and the coupling with 
the continuum becomes important. We shall focus here on the study of this coupling since, 
to our knowledge, the effect of the continuum coupling on thermodynamical properties of 
superfluid nuclei has not been investigated so far. 

The main problem in dealing with the continuum coupling in statistical calculations of 
excited nuclei is the fact that the particles moving in the continuum have a finite probability 
to be emitted from the nucleus. In other words, such processes are time dependent and they 
are difficult to accomodate in stationary models as BCS or HFB. But to extent stationary 
many body theories to time-dependent formalisms is not an easy undertaking. In fact it 
may even be a not well defined task since the initial conditions may induce chaotic solutions. 
These features were already recognized in the beginning of quantum mechanics. Thus, 
Gamow and Wigner tried to reconciliate the outgoing character of the decaying process with 
the conveniences of stationarity by solving the Schrodinger equation with outgoing boundary 
conditions (for references see ||). The corresponding solutions are related to the complex 
poles of the S-matrix, which define the so-called Gamow resonances. If the resonances are 
narrow the real parts of the complex poles of the S-matrix give the positions of the resonances 
while the corresponding imaginary parts give the decay widths. The narrow resonances are 
very important to describe nuclear correlations, especially in unbound or excited nuclei, 
because in these states the nucleons could move within the nuclear volume during a certain 
minimum time, so that they can interact with each other. One thus expects that a basis 
formed by bound and narrow Gamow resonances would provide a convenient framework to 



describe decaying processes |10| , |il| . However, the drawback of the calculations based on this 
representation is that one gets complex probabilities which are not always easy to interpret 

0- 

Another alternative to describe unstable nuclei is to use a basis consisting of scattering 
states, instead of Gamow states, in the vecinity of the resonant poles |12| , |13[ . In this case 
all quantities are real and one does not have problems of interpreting complex probabili- 
ties. Within this representation the widths of the resonances are obtained by evaluating the 
derivative of the corresponding phase shifts. This also defines the continuum level density 
rj| commonly used to estimate the contribution of the continuum to nuclear partition func- 



tions fl5H20||. The escape of the particles which move in resonant states is thus treated in 
these calculations as a stationary process, reflected by a constant particle density at large 
distances from the nucleus . The underling picture is an excited nucleus in dynam- 



ical equilibrium with an external nucleonic gas, whose contribution should be eventually 
extracted. Recently a similar framework was used to include continuum coupling in BCS 



equations at zero temperature |T2,13|. In this paper we will extend this method to study 



temperature dependent BCS equations (TBCS). 



II. FORMALISM 

The standard procedure to derive the temperature dependent BCS equations is to min- 
imize the grand potential corresponding to a pairing Hamiltonian. For a bound single 
particle spectrum with energies €j and a constant pairing interaction of strength G, the gap 
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and particle number equations at a finite temperature T are given by [|TJ 

2 _ v 1 - 2. 
G 2? 2E 3 



^ = E^[l-^(l-2/ J )] (2) 

where £y = [(e 3 - — A) 2 + A 2 ] 1 / 2 is the quasiparticle energy, fj = (l + exp(f3Ej))~ 1 is the Fermi 
distribution function and (3 = 1/kT. 

In principle the single particle energies e 3 - depend also upon temperature because the 
average mean field is a function of the nuclear excitations. However, in self-consistent 
temperature dependent Hartree-Fock calculations one finds |19| that for temperatures below 
T = 1 MeV, which is the range explored in temperature dependent BCS calculations, the 
single-particle spectrum is virtually the same as the one at zero temperature. We will also 
assume here that the single particle energies are those at zero temperature. 

The contribution of the continuum on thermodynamical properties of finite nuclei was 



studied mainly in connection with the problem of the liquid-gas phase transition [|Tq — |2C ] 
as well as in temperature dependent shell corrections [^TJ, but without including pairing 
correlations. In these calculations the effect of the continuum was introduced into the 
thermodynamical quantities through the level density. Thus the grand potential for a non- 
interacting system was taken as 



« = -T j Gfc(e) + g(e))ln(l + exp[-/3(e - \)])de (3) 
where <7&(e) = J2j <K e — e j) * s the level density of the bound spectrum and g(e) is the level 



density associated with the positive energy spectrum. In Ref. |19[ it is shown that the grand 
potential (|3]) describes a nucleus in dynamical equilibrium with a nucleonic gas. As discussed 
above, this is due to the fact that in a stationary treatment the nucleons scattered in the 
continuum are permanently emitted from the nucleus. To obtain the proper grand potential, 
i. e. the one corresponding to the nucleus itself, one should take away from Eq. @ the 
contribution of the nucleonic gas. This can be done by subtracting from the grand potential 
(0) the grand potential of the free nucleonic gas |19|], or by replacing in Eq. (|]) the level 
density g(e) by the quantity 



g(t) =9~ 9free = ~ ^ ^ (4) 

where gf ree is the level density in the absence of the mean field and 5j is the phase shift. The 
quantity g(e) is the continuum level density [|H |. It takes into account the contribution of the 



resonant part of the continuum spectrum. The continuum coupling can thus be included by 
replacing in the grand potential (|3|) the density g by the continuum level density g(e). This 
is a general recipe which can be applied to all quantities derived from the grand potential, 
as e. g. the energy and entropy of the system. 

The incorporation of the continuum in interacting systems can readly be performed fol- 
lowing a similar prescription. Thus the contribution of the continuum to the grand potential 
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of an excited superfluid nucleus can be expressed in term of the continuum level density as 
in Eq. (||), with the difference that now instead of the single particle energies one should 
use the quasiparticle energies. The corresponding TBCS equations can be obtained directly 
from Eqs. ([I]) and (fj) by replacing the level density of the bound states with the total level 
density, i. e. by g&(e) + g(e). Thus the TBCS equations with continuum coupling become, 

2 ^ l-2/ 3 - , f , , l-2/(e) J 

G = ^^E- + J 9{e) -2EW de > (5) 

N = E ^[1 - " 2/i)] + / g{e)\[l - ^(1 - 2f(e))]de (6) 

where the second term gives the contribution of the continuum to the pairing correlations. 
In the limit T=0 one gets the same equations as in Refs |12| , [r3| . For temperatures higher 
than the critical temperature the gap vanishes and the particle number equation is similar 



to the one used in thermodynamical calculations of non-interacting systems [|18| -|2T| . 

The contribution of the continuum to the energy and to the entropy can be introduced 
in a similar fashion, as mentioned above. One gets, 

A 2 



E = E n i e i + J g{e)n{e)ede - — , (7) 



S = EW/i) + MM} + I g(e){ln(Me)) + (3f u (e)E u (e)}de (8) 
j 

where rij is the occupancy of the state of energy €j, given by 

^[1-^(1-2/;)] (9) 

and fj = 1 + exp(—/3Ej). A similar notation is used for the corresponding quantities in the 
continuum. 

In the limit of vanishing widths the continuum level density becomes a sum of Dirac delta 
functions and the resonances act as bound states ("quasibound states"). This is the case for, 
e. g., protons trapped in a high Coulomb barrier. If a resonance is not narrow then a pair 
scattered to that resonance have a large probability to escape from the system. Therefore, 
its contribution to the pairing correlation is small as compared with the corresponding 
contribution from a quasibound state with similar energy and angular momentum |T2]JT3|| . 

Although narrow resonances play a fundamental role in the enhancement of pair cor- 
relations (and in most other measurable physical processes as well) their proper inclusion 
in the applications may be difficult. The reason for this is that in the region of a narrow 
resonance the level density increases abruptly and the numerical evaluation of the integrals 
in the equations above may require extremely small mesh intervals. One can circumvent this 
problem by noticing that the resonance is narrow because there is a pole of the S matrix 
which is very close to the real energy axis. Therefore, one can evaluate the integrals by 
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changing the integration path, that is by choosing a contour C in the complex energy plane 
which embodies the real axis around the narrow resonance, and by thereafter applying the 
Cauchy theorem Eq. (0) thus becomes 



r A 2 
E = ni€i + 53 n(E v )E v + g(e)n(e)ede - — 



(10) 



where n(E u ) are the occupation probabilities calculated in the complex poles E„ enclosed 
by the path C . If one neglects the contribution of the integral over the contour C in Eq. 
([10|) then the energy E would become complex. This is the case in representations based 
on Gamow resonances only [22] . As already mentioned, within such representations one 
obtains complex physical quantities which have to be interpreted. For the energy such a 
task is rather easy, since already by looking at the temporal evolution of the wave function 
one realizes that the real part corresponds to the actual energy of the system while the 
imaginary part is related to the corresponding decay probability. This interpretation is 
valid if the resonance is narrow, i. e. if the ratio between the width and the energy of 
the resonance is small p3[p4| . But the interpretation of complex probabilities in general is 
not so straightforward [n,Z5|. These problems do not appear here because in Eq. (ffH) the 
contribution of the integral over the contour C is included and therefore the energy E is 
real. 



III. NUMERICAL APPLICATION 

In order to illustrate how the continuum affects the properties of superfluid nuclei close 
to the drip line, in what follows we present the results given by the TBCS equations for the 
isotope 84 Ni. 

We calculated the single particle spectrum by using the HF approximation and the 



Skyrme III interaction [25|. The resonant energies are defined as the energies where the 



phase shift passes through tt/2 with a positive slope |2"3" |. The width is extracted from the 
value of the energy where the derivative of the phase shift is half of its maximum value. 
The bound states outside the closed shell N=50 and the resonant states considered in the 
TBCS calculations are listed in Table 1. It can be seen that these states form the equivalent 
of the major shell N=50-82. One notices that the relative positions of the single-particle 
states differ substantially with the ones corresponding to beta stable nuclei. Thus, the states 
with low angular momenta are shifted down as compared with the ones with high angular 
momenta. This is a general feature related to the diffusivity of the mean field in nuclei close 



to the drip line ||26|| , which is larger than the corresponding one in stable nuclei . 

The TBCS equations (Q) and (^) are solved starting with the HF spectrum calculated at 
zero temperature. In the absence of experimental information on heavy Ni isotopes in the 
open shell = 50 — 82 (the heaviest known Ni isotope is the double magic nucleus 78 Ni) 
we use for the strength of the pairing force the standard value G = 25 j A MeV, where A is 
the mass number. 

The variation of the gap with temperature is shown in Fig. 1. In the limit of T = 
the gap has the value A(0) = 0.955 MeV and decreases to zero at the critical temperature 
T c = 0.524 MeV. This critical temperature is smaller by about 3.8% than the one obtained 
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from the relation T c = 0.57A(0), which is the value predicted by a constant level density 
approximation 0. 

In Fig. 1 it is also shown the gap one would get if the width of the resonances are 
neglected, i.e. if in Eqs. (|)-|D one would take instead of the continuum level density a sum of 
Dirac delta functions. The gap at zero temperature is in this case bigger than in the previous 
case [ 12|,13| . As discussed above, the wider is a resonance the smaller is its contribution to 



the pairing correlation. That is, a pair scattered to a wide resonance spends less time in 
that state as compared to the time that the same pair would spend in a quasibound state 
and, therefore, the wider resonance contributes less to the pairing correlations. This has 
the important consequence that neglecting the widths of the resonances one increases the 
correlations of the system and therefore increases the critical temperature. This is what 
happens, e. g., with calculations that quantize the continuum by using an impenetrable 
box if the dimensions of the box are not extremely large (so that one gets a dense enough 
spectrum in the energy regions of resonant states). However, one has to mention that this 
is not a serious problem if the resonances are very narrow. 

Neglecting the widths one gets for the critical temperature the value T c = 0.722 MeV, 
which is only 1% larger than the value one would obtain in a constant level density approxi- 
mation. This indicates that this approximation would work quite well for nuclei close to the 
proton drip line, where the width of the proton single-particle resonances may be so narrow 
that they can be neglected. 

In Fig. 2 we show the dependence of the excitation energy upon temperature up to the 
point where the superffuid phase vanishes. It can be seen that if the effect of the width of 
resonant states is neglected, then the slope of the excitation energy becomes much smaller. 
This is also a manifestation of the fact that by neglecting the widths, the ground state 
becomes more correlated and therefore stiffer to thermal excitations. The same effect is 
observed for the entropy, as seen in Figure 3. 

At critical temperature the excitation energy has the value 2.185 MeV and the entropy 
is 5.708. For a constant level density g the excitation energy at critical temperature is given 
by m E c « aT c 2 + 0.5gA(0) 2 , where a = ^-g is the level density parameter, which in this case 
acquires the value a = 3.958 MeV -1 . We found that neglecting the widths the level density 
parameter increases to the value a = 4.143 MeV -1 . This implies that the difference in the 
excitation energies seen in Fig. 2 is esentially due to the mean occupancy of the resonant 
states (which is larger if the widths are neglected) and not due to the effective level density 
parameter. 

In conclusion, in this paper we have extended the temperature dependent BCS equations 
by introducing the coupling with the single-particle continuum. The contribution of the 
continuum is given by the resonant states and their effect is taken into account through the 
continuum level density. 

We found that the widths of the resonances affect significantly all physical quantities. In 
particular, the pairing correlations are diminished and this modifies significantly the value 
of the critical temperature at which the supefluid phase disappears. Also the dependence 
upon temperature of the excitation energy and of the entropy is considerably affected by 
the widths of the resonances, i.e. by their lifetime. 

However, in the case of proton superfluidity the single-particle resonant states close to 
the continuum threshold may have a very small width and therefore they can be treated in 
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the TBCS calculations as quasibound states, neglecting their widths altogether. 
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TABLES 



TABLE I. Neutron single particle states corresponding to Ni. These are the states used 
in the TBCS calculation. E n is the energy and L n the width of the state labelled by n. Both 
quantities are in MeV. 



n En L ra 

24/2 ~ -1-470 

3s 1/2 -0.730 

d 3/2 0.486 0.112 

g 7/2 1.605 0.010 

fen/2 3.296 0.016 



FIGURE CAPTIONS 

Figure 1: Dependence of the gap parameter A upon the temperature T. The dashed line 
corresponds to the case when the resonant states are considered as quasibound states. 

Figure 2: Excitation energy plotted as a function of temperature. The dashed line 
corresponds to the case when the resonant states are considered as quasibound states. The 
energy is plotted up to the critical temperature. 

Figure 3: Entropy plotted as a function of temperature. The dashed line corresponds 
to the case when the resonant states are considered as quasibound states. The entropy is 
plotted up to the critical temperature. 
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